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We compute spectra, waveforms, angular distribution and total gravitational energy of the gravi- 
es ; 

\f} . tiational radiation emitted during the radial infall of a massless particle into a Schwarzschild black 

hole. Our fully relativistic approach shows that (i) less than 50% of the total energy radiated to 
infinity is carried by quadrupole waves, (ii) the spectra is flat, and (iii) the zero frequency limit 
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agrees extremely well with a prediction by Smarr. This process may be looked at as the limiting 



case of collisions between massive particles traveling at nearly the speed of light, by identifying 



the energy E of the massless particle with mo7, ttiq being the mass of the test particle and 7 



the Lorentz boost parameter. We comment on the implications for the two black hole collision at 
nearly the speed of light process, where we obtain a 13.3% wave generation efficiency, and compare 



our results with previous results by D'Eath and Payne. 
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The study of gravitational wave emission by astrophysical objects has been for the last 
decades one of the most fascinating topics in General Relativity. This enthusiasm is of 
course partly due to the possibility of detecting gravitational waves by projects such as 
GEO600 ]T|], LIGO M or VIRGO ||, already operating. Since gravity couples very weakly 
to matter, one needs to have powerful sources of in order to hope for the detection of the 
gravitational waves. Of the candidate sources, black holes stand out naturally, as they 
provide huge warehouses of energy, a fraction of which may be converted into gravitational 
waves, by processes such as collisions between two black holes. As it often happens, the most 
interesting processes are the most difficult to handle, and events such as black hole-black 
hole collisions are no exception. An efficient description of such events requires the use of 
the full non linear Einstein's equations, which only begin to be manageable by numerical 
methods, and state-of-the-art computing. In recent years we have witnessed serious progress 
in this field [|J, and we are now able to evolve numerically the collision of two black holes, 
provided their initial separation is not much larger than a few Schwarzschild radius. At the 
same time these numerical results have been supplemented with results from first and second 
order perturbation theory [[J, which simultaneously served as guidance into the numerical 
codes. The agreement between the two methods is not only reassuring, but it is also in 
fact impressive that a linearization of Einstein's equations yield such good results (as Smarr 
puts it, "the agreement is so remarkable that something deep must be at work"). In 
connection with this kind of events, the use of perturbation methods goes back as far as 
1970, when Zerilli [f|J and Davis et al || first computed the gravitational energy radiated away 
during the infall from rest at infinity of a small test particle of mass mo, into a Schwarzschild 
black hole with mass M. Later, Ruffini || generalized these results to allow for an initial 
velocity of the test particle (this problem has recently been the subject of further study 
|1(J, in order to investigate the question of choosing appropriate initial data for black hole 
collisions). Soon after, one began to realize that the limit mo — > M describing the collision 
of two black holes did predict reasonable results, still within perturbation theory, thereby 
making perturbation theory an inexpensive tool to study important phenomena. 

In this paper we shall extend the results of Davis et al || and Ruffini |9j by considering a 
massless test particle falling in from infinity through a radial geodesic. This process describes 
the collision of an infalling test particle in the limit that the initial velocity goes to the speed 
of light, thereby extending the range of Ruffini's results into larger Lorentz boost parameters 
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7S. If then one continues to rely on the agreement between perturbation theory and the fully 
numerical outputs, these results presumably describe the collision of two black holes near 
the speed of light, these events have been extensively studied through matching techniques 
by D'Eath (a good review can be found in his book []TT]| ) , and have also been studied in 
|T2|]. The extension is straightforward, the mathematics involved are quite standard, but 



the process has never been studied. Again supposing that these results hold for the head 
on collision of two black holes travelling towards each other at the speed of light, we have 
a very simple and tractable problem which can serve as a guide and supplement the results 
obtained by Smarr and by D'Eath. Another strong motivation for this work comes from the 



possibility of black hole formation in TeV-scale gravity [13]. Previous estimates on how this 



process develops, in particular the final mass of the black hole formed by the collision of 



relativistic particles have relied heavily upon the computations of D'Eath and Payne [JIB 
A fresher look at the problem is therefore recommended, and a comparation between our 
results with results obtained years ago [11, 12) and with recent results [14, 15| are in order. 



Our fully relativistic results show an impressive agreement with results by Smarr [|1^] for 
collisions of massive particles near the speed of light, namely a flat spectrum, with a zero 
frequency limit (ZFL) (ttj)uj=o = 0.4244mQ7 2 . We also show that Smarr underestimated 
the total energy radiated to infinity, which we estimate to be AE = 0.26mQ7 2 /M, with M 
the black hole mass. The quadrupole part of the perturbation carries less than 50% of this 
energy. When applied to the head on collision of two black holes moving at the speed of 
light, we obtain an efficiency for gravitational wave generation of 13%, quite close to D'Eath 
and Payne's result of 16% fy], |16 |. 



Since the mathematical formalism for this problem has been thoroughly exploited over the 
years, we will just outline the procedure. Treating the massless particle as a perturbation, 
we write the metric functions for this spacetime, black hole + infalling particle, as 

where the metric g^b\ xV ) IS the background metric, (given by some known solution of Ein- 
stein's equations), which we now specialize to the Schwarzschild metric 

dv 2 

ds 2 = -f{r)dt 2 + J^ + r W + sin 2 6d<j) 2 ) , (2) 
where f(r) = 1 — 2M/r. Also, h a b(x u ) is a small perturbation, induced by the massless test 



particle, which is described by the stress energy tensor 

TV = -T^TyS J dX5 ^ x ~ Z W)^^ • (3) 

Here, z v is the trajectory of the particle along the world-line, parametrized by an affine 
parameter A (the proper time in the case of a massive particle), and po is the momentum 
of the particle. To proceed, we decompose Einstein's equations G a b = 8nT ab in tensorial 



spherical harmonics and specialize to the Regge- Wheeler |T7[ gauge. For our case, in which 
the particle falls straight in, only even parity perturbations survive. Finally, following Zer- 
illi's |L8| prescription, we arrive at a wavefunction (a function of the time t and radial r 



coordinates only) whose evolution can be followed by the wave equation 

d2Z ^ r) + y - V{r)\ Z(w, r) = (1 - 2M/r)S , (4) 

Here, the /—dependent potential V is given by 

_ f(r) [2(j 2 ((T + l)r 3 + 6a 2 r 2 M + 18arM 2 + 18M 3 ] 
V[r) ~ r^M + ar) 2 ' (5) 

where a = ^~ 1 ^ ;+2 - > anc l the tortoise coordinate r* is defined as = f(r). The passage 
from the time variable t to the frequency u has been achieved through a Fourier transform, 
Z(u;, r) = 73^172 f^lo e luJt Z(t, r)dt. The source S depends entirely on the stress energy tensor 
of the particle and on whether or not it is massive. The difference between massive and 
massless particles lies on the geodesies they follow. The radial geodesies for massive particles 
are: 

dT E dt E 



(6) 



dr f(r)(E 2 - 1 + 2M/r)V2 ' dr f(r) ' 
where E is a conserved energy parameter: For example, if the particle has velocity 
at infinity then E = ^_J, ^ 1/2 = 7. On the other hand, the radial geodesies for massless 
particles are described by 

dT 1 dt eo , x 

where again €q is a conserved energy parameter, which in relativistic units is simply p . We 
shall however keep eo for future use, to see more directly the connection between massless 
particles and massive ones traveling close to the speed of light. One can see that, on putting 
p — > m , e — > 7 and 7 — > 00, the radial null geodesies reduce to radial timelike geodesies, 



so that all the results we shall obtain in this paper can be carried over to the case of 
ultrarelativistic (massive) test particles falling into a Schwarzschild black hole. For massless 
particles the source term S is 

4^ e-^eo(4/ + 2) 1 / 2 a 

w(3M + ar) 2 ' [) 



To get the energy spectra, we use 

dE 1 (7 + 2)! .,- 

V|Z(cu,r)| 2 , (9) 



du 32tt(/-2)! 

and to reconstructe the wavefunction Z(t,r) one uses the inverse Fourier transform 

1 r°° 

Z(t,r) = j—tji J e**Z{u>,r)du>. (10) 

Now we have to find Z(u,r) from the differential equation @. This is accomplished by 
a Green's function technique. Imposing the usual boundary conditions, i.e., only ingoing 
waves at the horizon and outgoing waves at infinity, we get that, near infinity, 



Z = JL jf ZL Sdr. (11) 
Here, z L is a homogeneous solution of (§) which asymptotically behaves as 

z L ~ e"^* , r* -> -oo (12) 
2 L ~ B{u)e iur * + C(^)e-^ r * , n -> +oo. (13) 

W is the wronskian of the homogeneous solutions of @. These solutios are, zl which has 
just been defined, and zr which behaves as zr ~ e wr * , r* — > +oo. From this follows that 
W = 2iuC(uj). We find C(a;) by solving @ with the right hand side set to zero, and with 
the starting condition z^ — e~ tbjr * imposed at a large negative value of r*. For computational 
purposes good accuracy is hard to achieve with the form fll3|), so we used an asymptotic 
solution one order higher in \/{ur): 

+C(u)(l - ^l))e-^ ,r„ ^ +oo. (14) 

In the numerical work, we chose to adopt r as the independent variable, thereby avoiding the 
numerical inversion of r*(r). A fourth order Runge-Kutta routine started the integration of 




FIG. 1: Waveforms for the three lowest radiatable multipoles, for a massless particle falling from 
infinity into a Schwarzschild black hole. Here, the wavefunction Z is measured in units of eoPo- 

zl near the horizon, at r = r, = 2M + 2Me, with tipically e = 10~ 5 . It then integrated out to 
large values of r, where one matches Zl extracted numerically with the asymptotic solution 
([HP, in order to find C(u>). To find Z(t,r) the integral in fllOD is done by Simpson's rule. 
For both routines Richardson extrapolation is used. The results for the wavefunction Z(t, r) 
as a function of the retarded time u = t — r* are shown in Figure 1, for the three lowest 
radiatable multipoles, / = 2, 3 and 4. As expected from the work of Ruffini || and Ferrari 
and Ruffini , the wavefunction is not zero at very early times, reflecting the fact that the 



particle begins to fall with non zero velocity. At late times, the I = 2 (for example) signal 
is dominated by quasinormal ringing with frequency to ~ 0.35/M, the lowest quasinormal 
frequency for this spacetime [P0| . The energy spectra is shown in Figure 2, for the four 
lowest radiatable multipoles. First, as expected from Smarr's work ]T2 , the spectra is flat, 



up to a certain critical value of the frequency, after which it rapidly decreases to zero. This 
(/-dependent) critical frequency is well approximated, for each Z-pole, by the fundamental 
quasinormal frequency. In Table 1, we list the zero frequency limit (ZFL) for the first ten 
lowest radiatable multipoles. 

For high values of the angular quantum number /, a good fit to our numerical data is 

dE,\ 2.25 2 2 
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TABLE I: The zero frequency limit (ZFL) for the ten lowest radiatable multipoles. 
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ZFL (Xp-jr) 
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0.265 
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0.0068 


3 


0.075 


8 


0.0043 


4 


0.032 


9 


0.003 


5 


0.0166 


10 


0.0023 


6 


0.01 


11 


0.0017 



We therefore estimate the zero ZFL as 



( §Ei 
\ du) 



)uj=0 



W=ll (dE t \ 
1^1=2 \ dco Juj=0 



, 1 2.25 2^2 
T 2 12 2 toPo 



0.4244e^ . 



(16) 



To calculate the total energy radiated to infinity, we proceed as follows: as we said, the 
spectra goes as 2.25 /I 3 as long as uj < uj t Q N , where ujiq N is the lowest quasinormal frequency 
for that /-pole. For uj > u>iqn, dE/du ~ (In fact, our numerical data shows that dE/du ~ 
e -27awM ^ w ith a a factor of order unity, for uj > ujiqn)- Now, from the work of Ferrari 



and Mashhoon [^TJ and Schutz and Will |p2|| , one knows that for large /, ujiqn ~ 33/2 m • 
Therefore, for large I the energy radiated to infinity in each multipole is 

2.25(Z + 1/2) elpl 



A F 

1 " 3 3 / 2 / 3 M 
and an estimate to the total energy radiated is then 



(17) 



AE 



0.262 



foPo 
M 



(18) 



Let us now make the bridge between these results and previous results on collisions between 
massive particles at nearly the speed of light || [11], [12] . As we mentioned, putting p — > m 
and eo — > 7 does the trick. So for ultrarelativistic test particles with mass mo falling into a 

2 2 

Schwarzschild black hole, one should have (dE/du)^^ = 0.4244mg7 2 and AE = 0.262^^- 
Smarr 1T2T1 obtains 



M 



mtT = lrnh 2 - 0.4244mgr» 



Smarr 



0.2- 



M 



(19) 




FIG. 2: The energy spectra for the four lowest radiatable multipoles, for a massless particle falling 
from infinity into a Schwarzschild black hole. 
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FIG. 3: The energy radiated per solid angle as a function of 8. 

So Smarr's result for the ZFL is in excellent agreement with ours, while his result for the 
total energy radiated is seen to be an underestimate. As we know from the work of Davis 
et al P| for a particle falling from infinity with = most of the radiation is carried by 
the I = 2 mode. Not so here, in fact in our case less than 50% is carried in the quadrupole 
mode (we obtain AE l=2 = 0.1^, AE l=3 = 0.0447^ )• Tms is reflected in the angular 
distribution of the radiated energy (power per solid angle) 



dE _ AE (l~ 2)! 
dQ l (l + 2)\ 



(20) 



which we plot in Figure 3. Compare with figure 5 of ||12| . 

If, as Smarr, we continue to assume that something deep is at work, and that these results 
can be carried over to the case of two equal mass black holes flying towards each other at 
(close to) the speed of light, we obtain a wave generation efficiency of 13 %. This is in close 
agreement with results by D'Eath and Payne |Q, [Tj|, who obtain a 16% efficiency (Smarr's 
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results cannot be trusted in this regime, as shown by Payne [0]). Now, D'Eath and Payne's 
results were achieved by cutting an infinite series for the news function at the second term, 
so one has to take those results with some care. However, the agreement we find between 
ours and their results lead us to believe that once again perturbation theory has a much 
more wider realm of validity. To our knowledge, this is the first alternative to D'Eath and 
Payne's method of computing the energy release in such events. 

This work was partially funded by Fundagao para a Ciencia e Tecnologia (FCT) through 
project PESO/PRO/2000/4014. V.C. also acknowledges finantial support from FCT 
through PRAXIS XXI programme. J. P. S. L. thanks Observatorio Nacional do Rio de 
Janeiro for hospitality. 
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